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ABSTRACT 

We have discovered a 2.5 Mpc (projected) long filament of infrared-bright galaxies connecting two of the 
three ~ 5 x lO''* M© clusters making up the RCS 23 19+00 supercluster at z = 0.9. The filament is revealed in 
a deep Herschel Spectral and Photometric Imaging REceiver (SPIRE) map that shows 250-500 /im emission 
associated with a spectroscopically identified filament of galaxies spanning two X-ray bright cluster cores. We 
estimate that the total (8-1000 /^m) infrared luminosity of the filament is Lir ~ 5 x 10'^ Lq, which, if due to 
star formation alone, corresponds to a total SFR~ 900 M© yr"'. We are witnessing the scene of the build-up 
of a > 10'^ M© cluster of galaxies, seen prior to the merging of three massive components, each of which 
already contains a population of red, passive galaxies that formed at z > 2. The infrared filament demonstrates 
that significant stellar mass assembly is taking place in the moderate density, dynamically active circumcluster 
environments of the most massive clusters at high-redshift, and this activity is concomitant with the hierarchical 
build-up of large scale structure. 

Subject headings: Galaxies: clusters: individual (RCS 231953+0038.0, RCS 232002+0033.4, 
RCS 231948+0030.1) — Galaxies: high-redshift — Galaxies: starburst — Infrared: 
galaxies — Submillimeter: galaxies 



1. INTRODUCTION 

The average rate of galaxy growth, measured by the vol- 
ume averaged star formation rate (SFR), was a factor of ~10 
highe r at z ^ 1 than it is today (e.g., iHopkins & BeacomI 
|2006|) . when the environments destined to become the most 
massive (lO'^M©) clusters were still in the process of assem- 
bly. What effect does the growth of large scale structure at 
z ^ 1 have on the star formation histories of galaxies bound to, 
or being accreted onto, such environments? For over 30 years 
it has been known that the fraction o f blue star- forming galax- 
ies in clusters was higher in the past (iButcher & OemIerll978[ 
|1984|) . and more recent work has extended star formation sur- 
veys of distant clusters into the infrared (IR) regime (impor- 
tant for tracking the total SFR), reinforcing the view that there 
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has been s trong evol ution in t he total SFR of clusters since 
z -- 1 (e.g- lGeacheFa l. 2006; B ai et al.l2009l) . To some extent 
the evolution seen in cluster populations tracks the field, but 
it is important t o consider the environm ental context of this 
evolution (e.g., iGilbank & Balo gh 2008), since star-forming 
galaxies within high-z clusters are destined to evolve into the 
'red sequence' in the cores of the descendants of such en- 
viron ments today (e.g., iGilbank et al.ll2008^ iPoggianti et al.l 
l2008h . 

It is becoming clear that clusters exist at z '^ 1 that already 
contain an established population of massive galaxies that 
formed their stellar populations quickly at much higher red- 
shifts (z > 2) when the cluster was in a much earlier state of 
collapse (e.g.. lPapovich et ^l2010h . While this rapid forma- 
tion episode puts in place the massive tail of the cluster galaxy 
population, the dwarf end of the red sequence undergoes sig- 
nificant evolution via the continuous accretion of satellite 
galaxies, with the bulk of such assimilation in the form of 
dwarf galaxies at z = 0, st eepening the faint-end slope of the 
luminosity function (e.g., De Lucia et alj|2004t IGilbank et alj 
|2008a). What is the intermediate stage of red-sequence evo- 
lution at z ~ 1 and what role, if any, does the assembly of the 
clusters themselves play on their member galaxies? For exam- 
ple, is star formation triggered by, or at least associated with, 
the peripheral environments (outlying groups and filaments) 
feeding those massive clusters? 

In this Letter we discuss observ ations perforraed w ith the 
ESA Herschel Space Observatory ( Pilbratt et al.l 201 Oh Spec- 
tral an d Photometric Imaging REceiver (SPIRE; iGriffin et al.l 
12010'). We have discovered a remarkable filament of far- 
IR bright galaxies linking two of three massive clusters in 
close proximity at z = 0.9 - the RCS 2319+00 supercluster 
(IGilbank et al.l2008bh . Our study reveals a glimpse of the pro- 
cess of galaxy evolution that is coeval with the birth of one of 
the most massive structures in the Universe; a text-book ex- 
ample of structure assembly in the hierarchical paradigm. We 
assume cosmological parameters of JIa = 0.73, il^ = 0.27, 
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and //o = 7 1 km s"' Mpc"' dSpergel et al.ll2003h . 



2. RCS 2319+00 

RCS 231953+0038.0 (hereafter RCS 2319+00) was de- 
tec ted in the first Red- Sequence Cluster Survey (RCS- 
1; iGladders & Yed l2005h. Subsequent follow-up observa- 
tions with Chandra jHicks et alj 120081) and an extensive 
optical and near-IR spectroscopic campaign (iGilbank et alJ 
|2008b; A. Faloon et al., 2012 in preparation; R. Yan et 
al., 2012 in preparation) have revealed a remarkable super- 
cluster system (RCS 231953+0038.0, RCS 232002+0033.4, 
RCS 231948+0030.1) comprising three distinct X-ray lumi- 
nous cores with Lx ^ 3.6-7.6 x lO'^'^ergs"', separated by 
< 3 Mpc in the plane of the sky and ^ 10 Mpc along the line 
of sight (assuming a Hubble flow). X-ray, strong lensing and 
virial mass estimates all imply individual cluster masses of 
~ 5 X IO'^^Mq. The close proximity of the three components, 
taken with the fact that the extended X-ray emission and den- 
sity profiles of the red-sequence members appear to be aligned 
strongly suggests that the clusters are most likely in the early 
stages of a three- way merger that will r esult in a > 10'^ Mq 
cluster by z = 0.5 ( iGilbank et al.ll2008bh . 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. Herschel data 

RCS 2319+00 was observed with Herschel SPIRE using 
the 'Large Map' mode on 2009 December 17 as part of the 
Herschel Le nsing Survey (OD 217; OBSID 1342188181; 
lEgami et al.ll2010.) . The scan direction was set to scan angles 
A and B, with the length and height of the map set to 4', yield- 
ing a map 17' x 17' in extent, with a total observing time of 
1 .7 hr and on-source integration time of 0.6 hr ( 17 s per pixel). 
The data were reduced with the latest version of th e Hersche l 
interactive processing environment (HIPE v7.3.0; IOttll2010D . 
including the use of the sigma-kappa deglitcher which showed 
improvement over the default routine. The maps have default 
pixel scales of 6, 10, and 14" at 250, 350, and 500 ^Jm, re- 
specti vely. The HIPE SUSSEXtractor task (ISavage & Olived 
l2007h was used to produce match-filtered maps (raw maps 
convolved with a point spread function (PSF) of 18, 25, and 
36" at 250, 350, and 500 /xm, respectively) and signal-to- 
noise ratio (SNR) > 3 source catalogs. The 1 a rms instru- 
mental noise in the match-filtered maps is ~1.1, 1.1, and 
1 .7 mJy at 250, 350, and 500 fim, respectively. When quoting 
flux uncertainties, we add in quadrature to the instrumental 
noise a nominal confusion noise of 2±5.8, 6.3, and 6.8 mJy at 
250, 350, and 500 ^m, respectively (Nguyen et al. 2010). The 
weighted SPIRE map means are consistent with mJy. 

3.2. Spitzer data 

RCS-1 imaging and photometry catalogs from the 
24/L(m Multiband Imaging Photometer for Spitzer (MIPS; 
iRieke et al.1 12004 are presented in T Webb et al. (2012, 
in preparation). The MIPS coverage is complete down to 
~ 96 /xJy; for comparison, a local IR luminous galaxy, or 
LIRG, would be > 1 14 /iJy at the cluster redshift. Spitzer In- 
fraRed Array Camera (IRAC: Fazio et al...200 4) four-channel 
imaging was obtained through programs 30940 and 50720 
and reaches MJ^ + 1.5 at z = 0.9 and covers approximately the 
same area as the 24 /im observati ons. Image processing was 
performed using IRACPROC (iSchuster. M ai'engo & Patted 
12006) . software developed to wrap the existing Mosaick- 
ing and Point-source Extraction (MOPEX) pipeline in IRAC 



mode and add IRAC-specific improvements to outlier re- 
jection. Source detection and aperture photometry were 
perfo rmed using the Picture Processing Package (PPP; lYed 
Il99l . 

4. AN INFRARED-BRIGHT FILAMENT 

All three SPIRE maps show a region of apparent 'enhanced' 
far-IR emission over the background - in the form of a com- 
plex of blended and confused point sources - spanning the 'A' 
and 'B' components of RCS 2319+00, but not extending into 
the core region of component 'A' . For our analysis, we have 
defined a region encompassing the filament, described by the 
zone within 90" of a locus spanning the A-B cluster compo- 
nents (see Fig.[T]l. We discuss the significance of this emission 
below, but note that we cannot rule out the possibility that it 
is the result of an overlap of the outer regions of clusters A 
and B, rather than a true filament. However, we note that no 
similar structure is seen between cluster cores B and C which 
are similarly aligned on the sky, and the filament does not 
appear to run directly between the A and B cores, but rather 
traces an easterly arc (in both spectroscopy and SPIRE emis- 
sion), suggestive of a real and unique structure. In addition, 
the existence of the structure is supported by coincident over- 
densities of spectroscopic cluster members and MIPS 24 /im 
number counts. 

Based on extensive spectroscopy (> 2000 spectra) over 
the RCS 2319+00 structure with the VLT VIMOS and Mag- 
ellan IMACS instruments (A. Faloon et al., 2012 in prepa- 
ration; R. Yan et al., 2012 in preparation), A. Faloon et 
al. (2012, in preparation) have identified 14 24 fim spectro- 
scopically confirmed cluster members in the filament with 
0.8700 < z < 0.9247 and (43 ± 18) < S24 < (511 ±21)/iJy. 
This redshift range was chosen to encompass the entire su- 
percluster structure, based on the spectroscopic overdensity 
in Faloon et al. (2012, in preparation). We find that the fila- 
ment has a ~2 x higher confirmed spectroscopic member den- 
sity compared with regions in the surrounding 10 x 10 arcmin^ 
'field' area, with a significance of 2.5 a. While our errors 
above do not take into account spectroscopic completeness, 
we have attempted to minimize any bias in our comparison 
by limiting the analysis to 'field' regions with similar spec- 
troscopic sampling (ie., ^ 4 slits arcmin"^) as the filament re- 
gion and by also removing cluster members within 1 arcmin 
(0.5 Mpc) of the cluster cores. A more detailed analysis will 
be presented by A. Faloon et al. (2012, in preparation). There 
is also an apparent 2± 3 cr hint of an overdensity of 24 /im 
sources in the filament (N24 = (11.4 ±0.6) arcmin"^) versus 
the immediate background (N24 = (9.6 ± 0.2) arcmin"^). 

To get a sense of the significance of the apparent far-IR en- 
hancement within the filament as revealed by Herschel, we 
can simply sum the flux densities of the SNR> 3 cr 250 /im 
sources. However, note that point source extraction is inher- 
ently insensitive to blended and confused SPIRE sources and 
to extended emission, and thus our estimate will only yield 
a lower limit to the true level of far-IR emission if all of the 
emission seen is originating in the filament structure. To ac- 
count for the apparent diffuse emission in the filament region, 
we plot a pixel flux histogram of the filament versus the back- 
ground, showing a striking excess of far-IR emission in the 
filament at all flux levels (see Fig. |2]i. We now sum the 29 
250 /im sources in the filament region and subtract a conserva- 
tive area-normalized background estimate (based on the flux 
density of SNR> 3 sources outside of the filament region), 
yielding 5250 = 281 ± 33mJy. Due to the incompleteness in 
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Fig. 1 . — The Herschel filament in RCS 2319+00. The white outline indicates the filamentary region defined for the analysis, spanning components A-B, and is 
overlaid on the 250 /im SPIRE map (panels a and b) and on an RGB composite of the 250, 350 and 500 /jm maps smoothed to the lowest resolution {panels c and 
d). Smoothed X-ray emission contours (with the point-sources removed) are overplotted in yellow {panels a and h) or in black (panels c and d). Panel h) Blue 
'x' symbols indicate the positions of the 41 confirmed 0.8700 < z < 0.9247 cluster members, and are encircled in red if they are 24 /.tm-detected. Panel d) The 
red contours trace the surface density of galaxies selected by the RCS selection of 'red-sequence' galaxies at z ~ 0.9. As pointed out in Gilbank et al. (2008b;), 
note that the major axis of the red-sequence clusters tends to point towards its nearest neighbour, and also notice the striking offset from the red-sequence defined 
cluster cores 'B' and 'C to the peaks of the X-ray emission. 



our source extraction, unsurprisingly we see only a ~ 2 cr hint 
of an overdensity in the filament (N250 = ( 1 . 1 ± 0.2) arcmin"-) 
versus the background (N250 = (0.7 ±0.05)arcmin"-). 

We now present a more quantitative analysis and interpre- 
tation of the level of star-formation enhancement in the far- 
IR filament using our collection of extensive follow-up spec- 
ti-oscopy in RCS 2319+00. 

4.1. Total infrared emission and star formation rate in the 
filament 



We can estimate the contribution of far-IR emission in 
the filament by summing the far-IR flux densities of indi- 
vidual IR-bright cluster members and then scaling an appro- 
priate galaxy spectral energy distribution (SED). To identify 
IR-bright cluster members we work directly with the MIPS 
24 /xm imaging which has a much smaller PSF (FWHM~ 6") 
than the SPIRE data and therefore allows for reliable cross- 
referencing with the spectroscopy. We then determine the ap- 
propriate far- IR SED shape by performing a standard st acking 
analysis (e.g. iMarsden et al.ll2009t iPascale et al.ll2009h to de- 
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Fig. 2. — Area-normalized 250 /^m pixel flux histograms showing a strik- 
ing excess of far-IR emission in the filament versus the background (pixels 
outside the filament region). 

rive the average three-band SPIRE flux densities of the 24 fim- 
detected cluster members in the filament, thus yielding the av- 
erage SED from 24-500 fim. 

The SPIRE imaging is much shallower relative to the 
MIPS imaging, and is only sensitive to Ljr > IO'^Lq ac- 
tivity at the cluster redshift - and indeed only three of the 
24 /im cluster members appear to be 'detected' at 250 /im 
(see Fig. [Til. We thus proceed to stack the flux densities 
from the SPIRE maps at the location of the 14 MIPS cluster 
members (which have an average 24 /im flux of ~ 170 /iJy), 
yielding significant stacked flux densities of 11.0± I.6mJy 
(7.0cr), 10.7 ± 1.7 mJy (6.2a), and 8.5±1.9mJy (4.6a) at 
250, 350, and 500 /im, respectively. We note that stack- 
ing the cluster members on the SPIRE maps using the much 
more precise MIPS positions should help to minimize any 
blending and confusion contamination from the SPIRE emis- 
sion of non-cluster members. Fig. |3] shows the stacked 
SPIRE flux densities in relation to a set of spectral templates 
(jDale & Helou 2002) and the median //er; c/ieZ-Astrophysica l 
Terahertz Large Area Survey (H- ATLAS; lEales et al.llMIOl) 
galaxy template w ith 1 x 10" L© <Ldust < 3 x lO" L© from 
ISmith et al](l20I2l) . redshifted to z = 0.9 and normalized to the 
250 /im stacked flux density. The far-IR colors are consistent 
with the 'cool' templates from the Dale & Helou (2002) li- 
brary, with S(,i)/S\oo < 0.26 rest-frame colors. The H-ATLAS 
template is also a reasonable fit to the average far-IR emis- 
sion, yielding LiR = (1. 9 ±0.3) x lO'^L©, w hich is consistent 
with the LIRG-class. Using the iKennicutll ri998i) SFR con- 
version, SFR(MQyr-')=4.5 x I 0"^'*(Li r /erg s"' ) (whi ch as- 
sumes a starburst age <100Myr and a ISalpetenlI955l initial 
mass function), yields an average SFR of (30 ± 5) M© yr"' . 

At z = 0.9, the MIPS 24 /im band is dominated by a com- 
plex of broad polycyclic aromatic hydrocarbon (PAH) emis- 
sion. When assessing the total Lir, as defined by the 8- 
1000/tm integrated emission, normalizing template SEDs to 
24 /im emission alone is subject to a high degree of uncer- 
tainty, due to the potential variation (intrinsic and redshifted) 
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Fig. 3. — SEDs redshifted to z = 0.9 and normalized to the average 
(stacked) 250 ^m SPIRE flux at the 14 24 /im cluster members in the fila- 
ment. Three SED templates have been plotted for illustrative purposes, in- 
cluding the two extremes of the Dale & Helou (2002) library (in terms of the 
56o/5ioo color), as well as the iSmith et alj (i2012) H-ATLAS galaxy template 
with 1 X 10"L,7) < Ldust < 3 X lO" L©. The 'cool' color SED templates 
from'Dale & Heloul <20^ , as well as the H-ATLAS SED are well matched 
to the average far-IR emission of MlPS-selected emitters in the filament. At 
z = 0.9, the 24 /im band is tracing the rest-frame 13 fira light, which is dom- 
inated by a complex of broad RAH emission, thus we expect some variation 
in SED template matches in this band, as is seen here. 

of the PAH contribution in the 24 /im band. Similarly, our lack 
of spectral coverage in the 30-100 /im range means we have 
poor constraints on the contribution of a hot dust continuum 
component, which could boost the luminosity in the 24 /im 
band (e.g., Fig.[3]l. This is evidenced by the range of 24 /im 
fluxes predicted by the 250 /im normalized templates, which 
straddle the average 24 /tm emission. Therefore, we take the 
ISmith et aLl (120121) template as a conservative estimate of the 
form of die RCS2319H-00 members' SEDs, and use this to 
calculate the total Lir of the filament. 

Since our spectroscopic follow-up of the MIPS cluster 
members is incomplete, we can only place a conservative 
lower limit to the total far-IR emission in the filament by 
summing the 250 /im emission contributed by the 14 24 /im- 
detected cluster members, yielding ^250 = I54±22mJy. As- 
suming each galaxy has the same SED shape, we use the 
H-ATLAS template to determine a total Lir = (2.6 ± 0.4) x 
10'^ Lq, corresponding to a total SFR within the filament of 
(460±70)Mq yr"'. We stress that these estimates are strictly 
lower limits due to the spectroscopic incompleteness of the 
galaxies in the vicinity of the filament and the requirement of 
a 24/im detection. For example, if we include the 27 cluster 
members that were not detected at 24 /im in the stack, we find 
a factor of 40% more far-IR emission. 

Rather than attempt to coiTect for spectroscopic complete- 
ness we now use a different approach to estimate the total IR 
emission. Earlier we estimated the statistical total IR excess 
emission within the filament compared to the average over the 
map using only the SPIRE data. This should provide a more 
complete picture of the level of star-formation activity in the 
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filament as this estimate does not suffer from spectroscopic 
incompleteness. The total (background subtracted) signal is 
^250 = 28 1 ± 33 mJy within an area of 26.3 arcmin^, which cor- 
responds to a total LiR = (4.8 ±0.6) X lO'^L© and SFR within 
the filament of - 850 ± IOOMq yr"'. 

5. DISCUSSION 

How can we interpret these observations? First we note 
that star formation in the filament galaxies is unlikely to be 
impeded by ram-pressure stripping, given (a) th e relatively 
low d ensity of the intergalactic medium (IGM; iDave et all 
1200 Ih and relative galaxy velocities compared to the cores of 
the clusters (ram-pressure scales with the IC M density and 
the square of the galaxy velocity; iGunn & G ott 1972); and 
(b) the relatively short duration for the star formation events 
(lOOsMyr) compared to the time required for stripping in 
these environments (up to several Gyr, e.g., [McCarthy et al.l 
|2008|) . Furthermore, conditions within intermediate density 
(i.e., filament) environments might be suitable for tidally in- 
ducing activity via close interactions and mergers between 
gas-rich systems being accre ted onto such environments. For 
example. iGeach et al.l(l201 11) find an enhancement in the local 
fraction of (UV/IR selected) star-forming galaxies in a similar 
filamentary environment surrounding the CI 0016+16 super- 
cluster at z = 0.55. Is the star formation in the RCS 2319+00 
supercluster filament significantly enhanced compared to the 
field? 

We have established that the total SFR occurring within the 
filament, as traced by the 250 /im emission, is significantly 
larger than the local background. On its own, however, this 
does not directly indicate triggered or enhanced star forma- 
tion as this may simply reflect the underlying over-density of 
star-forming galaxies which have not been actively quenched. 
To quantify any possible enhancement we look at the SFR per 
unit steflar mass, M*, or specific SFRs (SSFRs) of the MIPS 
cluster members in the filament. The stellar mass was derived 
directly from the measured IRAC-4.5 /im flux densities after 
app lying an (Rc-z') co\or correction to the M/L ratio, follow- 
ing iBeiretal] (l2003h . Normalized to the same stellar mass, 
the individual IR galaxies within the filament have ^3x the 
SSFR of IR-bright galaxies within the central cluster cores, 
but only ^ 1 .5 x that of galaxies in the rest of the supercluster 
To avoid systematic differences we do not attempt to compare 



this to field galaxies at this epoch and simply note that the fil- 
ament galaxies do not show enhanced emission compared to 
the rest of the supercluster 

How important would this star formation be to the overall 
growth of what will be a very massive cluster? Given the ap- 
parent low stellar mass of the individual galaxies within the 
filament, their eventual fate is to fall onto the faint-end of the 
red sequence. We have estimated the total stellar mass al- 
ready in place in RCS 23 19+00-A within 2. 1 Mpc to be M,^ -- 
3 X 10'^ M0. This estimate is obtained by fitting a Schechter 
function to the background-subtracted 4.5 /im flux and inte- 
grating this to provide a total rest-frame /T-band luminosity 
(above Mab ^ -22) which is then converted to stel lar mass 
assuming an early-type SED ( jLeitherer et al.|[l999l) . Using 
the Mtot derived from the X-ray luminosity provides a stellar- 
to-total mass conversion which we then apply to the other two 
cluster cores, yielding a total M*(<2.1 Mpc)-- 7 x lO'^M©. 
Thus, if the filament galaxies continue at their current SFRs 
until the clusters finally merge at z ^0.5, they will increase the 
descendent cluster mass by 20%. This should be regarded as 
an upper-limit, as the star formation will begin to shut down 
at earlier times as the galaxies enter the increasingly dense 
cluster environment and as fuel is exhausted; nevertheless the 
pre-processing witnessed in this filament appears to be an im- 
portant stage of the build-up of the cluster stellar mass. 
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